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Strongly interacting fermions define the properties of complex matter throughout nature, from 
atomic nuclei and modern solid state materials to neutron stars. Ultracold atomic Fermi gases have 
emerged as a pristine platform for the study of many-fermion systems. Here we realize a quantum gas 
microscope for fermionic 40 K atoms trapped in an optical lattice, which allows one to probe strongly 
correlated fermions at the single atom level. We combine 3D Raman sideband cooling with high- 
resolution optics to simultaneously cool and image individual atoms with single lattice site resolution 
at a detection fidelity above 95%. The imaging process leaves the atoms predominantly in the 3D 
motional ground state of their respective lattice sites, inviting the implementation of a Maxwell’s 
demon to assemble low-entropy many-body states. Single-site resolved imaging of fermions enables 
the direct observation of magnetic order, time resolved measurements of the spread of particle 
correlations, and the detection of many-fermion entanglement. 


The collective behavior of fermionic particles gov¬ 
erns the structure of the elements, the workings of 
high-temperature superconductors and colossal magneto¬ 
resistance materials, and the properties of nuclear mat¬ 
ter. Yet our understanding of strongly interacting Fermi 
systems is limited, due in part to the antisymmetry re¬ 
quirement on the many-fermion wavefunction and the 
resulting “fermion sign problem”. In recent years, ul¬ 
tracold atomic quantum gases have enabled quantita¬ 
tive experimental tests of theories of strongly interact¬ 
ing fermions urtu In particular, fermions trapped in 
optical lattices can directly simulate the physics of elec¬ 
trons in a crystalline solid, shedding light on novel phys¬ 
ical phenomena in materials with strong electron corre¬ 
lations. A major effort is devoted to the realization of 
the Fermi-Hubbard model at low entropies, believed to 
capture the essential aspects of high-T c superconductiv¬ 
ity [5HTT]. For bosonic atoms, a new set of experimen¬ 
tal probes ideally suited for the observation of magnetic 
order and correlations has become available with the ad¬ 
vent of quantum gas microscopes mm, enabling high- 
resolution imaging of Hubbard-type lattice systems at the 
single atom level. They allowed the direct observation 
of spatial structures and ordering in the Bose-Hubbard 
model H3H5] and of the intricate correlations and dy¬ 
namics in these systems mm- a longstanding goal 
has been to realize such a quantum gas microscope for 
fermionic atoms. This would enable the direct probing 
and control at the single lattice site level of strongly cor¬ 
related fermion systems, in particular the Fermi-Hubbard 
model, in regimes that cannot be described by current 
theories. These prospects have sparked significant experi¬ 
mental effort to realize site-resolved, high-fidelity imaging 
of ultracold fermions, but this goal has so far remained 
elusive. 

In the present work, we realize a quantum gas micro¬ 
scope for fermionic 40 K atoms by combining 3D Raman 
sideband cooling with a high resolution imaging system. 


The imaging setup incorporates a hemispherical solid im¬ 
mersion lens optically contacted to the vacuum window 
(Fig. 1(a)). In combination with a microscope objec¬ 
tive with numerical aperture (NA) of 0.60, the system 
achieves an enhanced NA of 0.87 while eliminating aber¬ 
rations that would arise from a planar vacuum window. 
In order to keep the atoms localized while performing flu¬ 
orescence imaging, one must simultaneously cool them in 
order to mitigate the heating from spontaneously emitted 
imaging photons. Previous microscope experiments in 
Hubbard-type lattices mm cool via optical molasses. 
In contrast, we employ 3D Raman sideband cooling he- 
[26], in which Raman transitions on vibration-lowering 
sidebands are combined with optical pumping to provide 
cooling. Our method therefore not only achieves site- 
resolved imaging, but also leaves a large fraction of the 
atoms (72(3)%) in the 3D motional ground state of each 
lattice site. This opens up prospects for the prepara¬ 
tion of low entropy many-body states, by measuring the 
atoms’ initial positions and rearranging them into the 
desired configuration m- 

Raman sideband cooling has previously been used to 
cool 87 Rb and 133 Cs atoms in lattices and in optical 
tweezers to large ground state populations [T9H26] . Here, 
we realize continuous Raman sideband cooling of 40 K 
using two states from the ground hyperfine manifolds, 

|a) = \F = 9/2,= -9/2) and | b) = |7/2, —7/2), 
which form an approximate two-level system. To make 
|a) and | h) non-degenerate with other hyperfine states, 
and to provide a quantization axis for optical pump¬ 
ing, we apply a bias field of 4.2 G along the x direction 
(Fig. 1(b)). A pair of Raman beams collinear with x and 
y lattice beams, but not retro-reflected, drives vibration¬ 
lowering Raman transitions from | a) to \b) (Fig. 1(c)). 
The Raman lasers are detuned —41 GHz from the D2 
transition. The optical pumping is performed on the D1 
transition, 3nm away from the D2 line, allowing us to 
filter out stray Raman light while transmitting atomic 
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|9/2, —9/2) 

FIG. 1. (a) High-resolution imaging setup. A solid immersion 
lens is formed by a spherical cap and a super-polished sub¬ 
strate contacted on either side of the vacuum window. Using 
an objective with NA = 0.60, the system achieves an effective 
NA s= 0.87. The substrate reflects 1064 nm light while trans¬ 
mitting D1 and D2 light of 40 K. The lattice beams are shown 
in red; the optical pumping and x-Raman beams are shown in 
blue, (b) Top view of Raman beams (blue) and optical pump¬ 
ing beam (green), (c) Raman cooling scheme. The Raman 
beams detuned near the D2 line (solid blue) drive vibration¬ 
lowering transitions. The optical T-pumping beam (dashed 
green) is tuned to the D1 line. Not shown is the raF-pumping 
beam. 


fluorescence. By collecting the photons that are spon¬ 
taneously scattered during this optical pumping process, 
we can image the atoms without using additional reso¬ 
nant light. 

To prepare a cold cloud of fermionic atoms under the 
microscope, 40 K is first sympathetically cooled with 23 Na 
in a plugged magnetic quadrupole trap [28] . centered 
^ 9 mm below a super-polished substrate that forms the 
bottom of the solid immersion lens. After removal of 
23 Na, the cloud of ^ 1 million 40 K atoms is magnetically 
transported to the substrate and trapped in a vertical lat¬ 
tice formed by a 1064 nm laser beam reflected off the sub¬ 
strate at an angle of 5.9°. A single layer 7.8 fi m from the 
surface is selected using a radiofrequency sweep in a verti¬ 
cal magnetic gradient followed by a resonant light pulse 
that removes atoms in the remaining layers. Next, we 
prepare a 50/50 mixture of |9/2, —9/2) and |9/2, — 7/2) 
to allow thermalization, and transfer the atoms to a ver¬ 
tical (z direction) 1064 nm beam, forming a lattice along 
z with a spacing of 532 nm (Fig. 1(a)). After evaporat¬ 
ing by lowering the power of the ^-lattice, the z-depth is 



FIG. 2. (a) Electronic ground (| g,N)) and excited (|e, TV)) 
states dressed by N optical pumping photons, (b) On res¬ 
onance, all dressed states are equally populated and expe¬ 
rience an anti-trapping potential. The four decay channels 
between the dressed states are equal, (c) The dressed states 
with Q/S = 0.175, shown with the dark state | D). Only a 
small fraction fl 4 /(5 4 of the steady-state population resides in 
anti-trapping states. For the atoms in the trapping state, the 
branching ratios between trapping transitions (solid arrows) 
and anti-trapping transitions (dashed arrow) are shown. 


increased to 180 /rK. We simultaneously ramp up two ad¬ 
ditional 1064 nm beams (Fig. 1(a)) reflected off the sub¬ 
strate at 10.8° and retro-reflected. These form a lattice 
in the horizontal plane with a spacing of 541 nm. 

During imaging, the atoms are trapped in a deep 
lattice, where the potential at each lattice site can be 
approximated by a harmonic well with vibrational fre¬ 
quency uj. At our imaging depth, the vibrational fre¬ 
quencies for the three axes are (w x ,u) y ,u) z ) = 27r x 
(280,300,260) kHz, corresponding to lattice depths of 
210/iK, 240 /iK and 180/iK respectively. The Rabi cou¬ 
pling for transitions that change the vibrational num¬ 
ber by one is proportional to the Lamb-Dicke param¬ 
eter, r] = Aka , where a = yis the harmonic 
oscillator length and hAk is the momentum transfer 
due to the Raman beams. Along our lattice directions, 
Ak x = 8.0 yum -1 , Ak y = 8.0 yum -1 and Ak z = 3.1 yum -1 , 
yielding Lamb-Dicke parameters of 0.17 for x and y, and 
0.068 for z. The polarizations of both Raman beams are 
linear and parallel to the substrate (Fig. 1(b)), in or¬ 
der to avoid differential effective magnetic fields between 
|a) and | b) that would arise for circularly polarized light. 
The Raman beam along the y -axis contains a single fre¬ 
quency, whereas the Raman beam along the x-axis con¬ 
tains three frequencies, allowing us to address the cooling 
sidebands of the three directions simultaneously. The re¬ 
sulting two-photon detunings from the bare | a) —» \b) 
transition are 400 kHz, 450 kHz and 360 kHz for cooling 
along x, y and 2 respectively. These frequencies compen¬ 
sate for differential Stark shifts that arise in the presence 
of optical pumping light. The r-Raman beam intensi- 
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FIG. 3. (a) Site-resolved imaging of fermonic atoms on a densely-filled 541 nm-period optical lattice, with an exposure time of 
two seconds; one can clearly discern the lattice structure and individual atoms, (b) Raman spectrum after cooling. The lattice 
depths are chosen such that the vibrational sidebands are well-resolved. The heating sidebands for the three lattice axes are 
labeled Z, X and Y. We observe a large sideband asymmetry, from which we extract a 3D ground state occupation of 72(3)%. 


ties of the three frequency components are 0.79 W/cm 2 , 
0.47W/cm 2 and 0.49 W/cm 2 respectively; the intensity 
of the y -Raman beam is 2.0 W/cm 2 . 

In addition to these Raman beams, optical pumping 
light is present to complete the cooling cycle. During op¬ 
tical pumping, atoms enter electronically excited states, 
and preferentially decay into the desired state. Typi¬ 
cally, the excited states experience an anti-trapping po¬ 
tential when the ground state experiences a trapping po¬ 
tential. For our 1064 nm lattice, the anti-trapping poten¬ 
tial for the atoms in the 4 P 1/2 states is 5.4 times stronger 
than the trapping potential for atoms in the 4 Si/ 2 states, 
due to the 4 P;l/ 2 3D 3 / 2 transition at 1169 nm. This 

strong anti-trapping would lead to heating and diffusion 
of atoms through the lattice during imaging. 

A solution to this problem is to detune the optical 
pumping light away from resonance. This reduces the 
population in the anti-trapping states and favors tran¬ 
sitions into trapping states. One can understand this 
by considering the dressed states of a driven two-level 
system. At zero intensity, the dressed states, | g,N) and 
|e, TV), correspond to the bare trapping and anti-trapping 
states with N photons respectively (Fig. 2(a)). In the 
presence of resonant pumping light, neither of the dressed 
states is trapping (Fig. 2(b)), which leads to heating. 

However, at large detunings S 12, where 12 is the 
Rabi frequency, one dressed state becomes trapping. Fur¬ 
thermore, spontaneous decay among the dressed states 
favors population in the trapping states. Specifically, the 
decay rates Y eg , r ee , Y gg and T ge , defined in Fig. 2(b), 
are proportional to 1, 8, s and s 2 respectively, where 
s = fl 2 /S 2 . The ratio of the anti-trapped population to 
the trapped population is suppressed, because in steady 
state, it is given by Y ge /Y eg which scales as s 2 . An¬ 


other benefit of large detunings involves the state | D) 
into which atoms are optically pumped (Fig. 2(c)). This 
state is dark to the optical pumping light, and hence 
has no excited state admixture; consequently, it expe¬ 
riences a trapping potential. Atoms in trapping states 
\g,N) decay preferentially into \D), since the ratio of 
anti-trapping transitions to dark state transitions scales 
as T ge /T g g = s, which is small at large (5. 

In light of these considerations, we detune the 
hyperfine-changing (F) pumping light —80 MHz from the 
Stark-shifted F = 7/2^F' = 9/2 transition and the 
Zeeman-level (mp) pumping light —80 MHz from the 
Stark-shifted F = 9/2 -T F' = 7/2 transition. The 
optical pumping beam co-propagates with the r-Raman 
beam, and has its polarization optimized for minimal cr + 
admixture. The F and mp components have intensities 
of 5.8 mW/cm 2 and 1.6 mW/cm 2 respectively; the Lamb- 
Dicke parameters for optical pumping are 0.17 for x and 
2/, and 0.18 for z. 

This Raman cooling scheme allows us to collect fluo¬ 
rescence while keeping the atoms confined to their lattice 
sites (Fig. 3(a)). Furthermore, we find that atoms are 
cooled predominantly into their motional ground state. 
Indeed, Raman spectroscopy reveals a 3D ground state 
population of 72(3)% after cooling (Fig. 3(b)). Note that 
parameters are optimized for imaging fidelity rather than 
for a large ground state population. We measure a fluo¬ 
rescence rate of ~ 5000 photons/atom • s with a lifetime 
of ~ 30 seconds. With a photon collection and detection 
efficiency of 20%, about 1000 photons per atom can be 
collected in an exposure time of one second, which is suf¬ 
ficient to detect single atoms with high fidelity. Using the 
number of scattered photons and the Lamb-Dicke param¬ 
eters for spontaneous emission, we estimate the cooling 
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FIG. 4. (a) Radially averaged PSF extracted from isolated atoms; the FWHM is 640 nm. (b) Intensity histogram after binning 
the deconvoluted images by lattice site. The threshold for reconstruction is shown by the dashed line, (c) Loss and hopping 
rates, shown in red squares and blue circles respectively, as extracted from 20 consecutive one-second exposures, (d) Correlation 
measurement g^ 2 \r) of a thermal cloud, showing the absence of distance-dependent loss, (e) Image of sparsely filled lattice with 
grid lines showing lattice spacing and orientation. The exposure time is one second. (f,g) The same image after deconvolution 
and with the filled sites identified. 


rate to be ~ 4 /iK/ms. 

To verify that we can resolve individual lattice sites, 
we measure the point spread function (PSF) of our imag¬ 
ing system using isolated atoms from sparsely populated 
images. The measured PSF has a full-width-half-max 
(FWHM) of 640 nm (Fig. 4(a)). Images are deconvolved 
with the PSF to achieve sub-lattice-site resolution. From 
such images, we also extract the lattice axes and spacings 
necessary to reconstruct the atomic distribution. Binning 
the intensity of the deconvolved image by lattice site re¬ 
veals a clear bimodal distribution (Fig. 4(b)), which is 
used to determine whether a site is filled (Fig. 4(e,f,g)). 
This bimodality gives a reconstruction error of < %1. 

An important aspect of quantum gas microscopy is the 
fidelity of the imaging process, which can be character¬ 
ized by hopping and loss rates. To this end, we take a 
series of images of the same atomic cloud and observe 
changes in the site occupations between images. Sites 
that are empty but become occupied in a subsequent im¬ 
age are counted as hopping events; sites that become 
empty are counted as loss events. The Raman cooling 


parameters are optimized for low hopping and loss rates 
while maintaining a fixed level of fluorescence. For opti¬ 
mized parameters, we achieve loss rates of < 4.8 db 0.2% 
and hopping rates of < 1.2 ± 0.2% for one-second ex¬ 
posures of clouds with fillings between 0.10 and 0.20 
(Fig. 4(c)). These rates, which include reconstruction 
errors, give an detection fidelity of > 95% for sparse 
clouds. At higher fillings, hopping events lead to loss of 
additional atoms in doubly-occupied sites due to light- 
assisted collisions. However, even for unity filling, we 
estimate the imaging fidelity to still be > 94%, because 
of the low hopping rate. 

To ensure that the imaging does not cause additional 
losses for neighboring atoms, one can measure the 2- 
point correlation function g^ 2 \r) = (n(x)n(x + r)) / ( n ) 2 
of thermal clouds, since distance-dependent loss will pro¬ 
duce anti-correlations at short distances. For a dilute 
thermal cloud with a filling of 0.19, one 29 x 32 site im¬ 
age gives g( 2 \r) = 1.00(7) for distances from r m 1 to 
10 lattice spacings, indicating that the imaging does not 
cause significant distance-dependent loss (Fig. 4(d)). 
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In conclusion, we have realized high-fidelity site- 
resolved imaging of 40 K fermionic atoms in a Hubbard- 
type optical lattice by combining 3D Raman sideband 
cooling with high resolution photon collection. In con¬ 
trast to existing boson microscopes, the technique leaves 
atoms predominantly in the absolute 3D ground state 
of a given lattice site. This opens up new ways to 
assemble low-entropy Fermi-Hubbard systems atom by 
atom [271 Ell [30]. Combining site-resolved imaging with 
on-site manipulation would allow one to deterministically 
create localized excitations and follow their time evolu¬ 
tion m Finally, the presence of 23 Na in our system 
invites the realization of a quantum gas microscope for 
ultracold fermionic NaK molecules Em. which have been 
proposed as a new resource for quantum information pro¬ 
cessing and quantum simulation of lattice models with 
long-range dipolar interactions. 
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